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We study steady-state spherically symmetric accretion of a galileon field onto a Schwarz- 
schild black hole in the test fluid approximation. The galileon is assumed to undergo a stage 
of cosmological evolution, thus setting a non-trivial boundary condition at spatial infinity. 
The critical flow is found for some parameters of the theory. There is a range of parameters 
when the critical flow exists, but the solution is unstable. It is also shown that for a certain 
^ range of parameters the critical flow solution does not exist. Depending on the model the 

sound horizon of the flow can be either outside or inside of the Schwarzschild horizon. The 
latter property may make it problematic to embed the galileon theory in the standard black 
hole thermodynamics. 

I. INTRODUCTION 
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Galileon is a scalar field theory, introduced in [1]. This theory has a remarkable property: 
its Lagrangian contains higher-order derivative terms, while the equations of motion are of only 
second derivatives of the field. The galileon theory is invariant under the galilean transformation 
of the field. It was originally introduced as a generalization of the scalar field theory, left over 
in the decoupling limit [2] of Dvali-Gabadadze-Porrati (DGP) model 0] and which allows the 
implementation of the Vainshtein mechanism jj]. The covariant version of galileon (covariant 
£NJ , galileon) was found in 0], where a model yielding only up to second-order derivatives in both field 

and metric was constructed. Other aspects of the galileon and related ideas were developed later 

a. 

The galileon and related models were also applied for construction of cosmological models [Tj 
(a more generic scalar-tensor model was suggested in 0| where the gravity is modified in infra-red, 
while the nonlinear kinetic coupling provides the restoration of general relativity). More recently, 
several galileon-related models were studied in the cosmological context, see e.g. 
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The galileon model is interesting in several aspects. First of all, this is a theory with non- 
quadratic kinetic coupling. This leads to propagation of perturbations in an effective metric, 
which is generically different from the gravitational one. In particular, the propagations may have 



sub- or superluminal velocities on non-trivial backgrounds. And indeed, as it was shown in [1|], 
galileon generically contains superluminal excitations (the same happens in for a scalar left over 
in the decoupling limit of DGP [3])- The existence of faster-than-light propagation on non-trivial 
backgrounds does not necessary im ply inconsistencies 14]. In particular, "superluminal" DBI 



inflation [15] or accreting k-essence [16[ are free of any time paradoxes (Note, that in [12], however, 
it was argued that superluminality of certain theories might indicate the problem for embedding 
them into a UV complete theory) . Another interesting feature of galileon was recently discussed in 
0, cosmological models with the phantom behavior were constructed, and it was argued that 
perturbations are ghost-free for some choice of parameters and initial conditions. 

In the original non-covariant version of the galileon it was found that the coefficient of the 
Lagrangian must be constrained, depending on the background [l|. The check was performed by 
constructing static spherically symmetric solutions with a source in the de Sitter universe. It is thus 
seems important to continue this route and to test the galileon model further against pathologies. 
In particular, the covariant version of the galileon must be checked for existence of solutions and 
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their stability. One of the useful area to test theories for unusual/pathological behavior is black 
hole physics. Indeed, previously the study of models in the black hole background has revealed 
interesting (and sometimes undesired) features of und erly ing models, e.g. it was shown that the 



phantom field accretion decreases the black hole mass [13] ; the ghost condensate plus a black hole 
forms a perpetuum mobile of the 2nd kind [3] (however in it was argued that ghost condensate 
does not violate the second law of thermodynamics); and the superluminal k-essence allows one to 
send signals from inside of the event horizon of a black hole [16]. 

In this paper we consider the galileon accretion onto a Schwarzschild black hole. The paper is 
organized as follows. In Sec. 2 we specify the action of the theory, give equations of motion and 
the energy-momentum tensor for different types of galileon. Section 3 is devoted to the galileon of 
the certain type, namely, containing only the terms left over in the decoupling limit of DGP. We 
construct the family of solutions describing a steady-state inflow of the scalar field. The unique 
physical solution is fixed by the requirement that it is regular both at the Schwarzschild and sound 
horizons. In Sec. 4 we widen our study to include another form of the galileon term. In Sec. 5 
we briefly summarize our results, discuss applicability of our study and physical consequences. 
In particular, for the galileon of the DGP type we find no constrains on the Lagrangian. When 
another galileon term is included, we obtain an additional constraint as compared to pj]. For a 
certain choice of parameters, the galileon forms the sound horizon inside the Schwarzschild one, 
thus giving an opportunity to send signals from the inside of the Schwarzschild horizon. We also 
discuss accretion of the "phantom" type of galileon. 

II. MODEL 

The general scalar field covariant galileon action reads [H] 

Sir = I d^xy^gL^, (1) 



where Lagrangian density for the can be written as a linear combination of terms, 

i=5 

Cji = ^ ^ Cj£j , (2) 
i=l 

where, 

£i = 7T, £ 2 = 7r ;At 7r ;M , £ 3 = 7r ;A1 7r ;M D7r, 
A = (7r ;a 7n 



2 (Ott) 2 - 2 (tt^ttH - \ V ^R 



(3) 



and the term £5 has a more complicated structure and contains higher order derivatives in ir, we 
are not going to study it in this papei0. We will also set c\ to zero in the following study, i.e. the 
"potential" term is not included. Otherwise, if c\ 7^ 0, there is no steady-state solution for the 
accretion (this is similar to the case of the k-esence field see, e.g. [ioj]). The constant c 2 in (|2j) is 
dimensionless, C3 has dimension —3, in the context of the DGP model this coefficient is usually 
written as r c /Mpj where r c is the so called cross-over scale and Mpi is the Planck mass. The 
coefficient C4 has mass dimension —6. 



1 Throughout this paper we use the signature (—,+,+,+). We also follow here the notations of Ref. [f|, which 
are different to ones in [J. Namely, neglecting gravity, and up to the integration by parts, £2 = — 2£ 2 VHT , 
£3 = —2C^ RT and £4 = — 4£ 4 vflT , where £f flT are the definitions of To avoid confusions with signs we 
also assume that our full action, including gravity, has the overall minus sign compared to so that, e.g. the 
constraint c?2 > in 1], for asymptotically flat space-time implies C2 > in our notations. 
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The energy- momentum tensor derived from J2]), ([3]) as T$ = 2/y/— g (SS^/dg^) is 

Z#) = 2 V> „ - g pu (dvr) 2 , (4a) 
^ )V Utt - 27r i(M V„) (0tt) 2 + ^7T' Q Va (5vr) 2 , (4b) 
T$ = -8 (Dtt) vr" [vr ;M ir. pv + vr ;!/ tt ;pp ] + 4 (Dvr) 2 (vr ;p 7r.„) - 4 (Dtt) (tt ; a ^ ;A ) (tt.^) 

- 8 (vr ;A tt' X p tt ;p ) (tt.^) + 8 (vn A ^ ;Ap ) (tt* ti>) - 4 (tt ; a p tp*") (vr ;p 7T ; „) 
+ 4 (vr ;A vr ;A ) (vr ;pp tt^) + 8 tt ;A tt*' [vr 

+ 8 (Dvr) (n. x n' Xp tt ;p ) - 8 (tt.* tt ]Xp n, pa vr ;,J ) ^ 

- 2 (vr ;A vn A ) (7r ;p(7 th^) g pv - 2 (tt ;A ^ ;A ) (vr ;p V ) J2 + i (vr ;A 7T A ) (tt ;p tt") 
+ 4 (vr ;A vr A ) vr^ [R pfl tt, v + ii pi , tt ;p ] - (tt ;A ^ ;A ) (vr ;p tp') 



- 4 (^ ;A 7T A J (vr ;p vr ;CT ) ^ + 4 (^ ;A 7P A J (vr" ir' a fl p/W(T ) . (4c) 

Note that the expression for Tm} differs from the one derived in [3] by factor —2, because of the 
difference in the definition for the energy-momentum tensoiH. Equations of motion are obtained 
by the variation of (pQ) with respect to tt, £u\ = (—g)~ l l 2 8Su\/5iT = 0, with 

£ 2 = -2D7T, 

£ 3 = -2 ((&) 2 - (VW) 2 - R^ VtV ) , 

£ 4 = -4 (Dtt) 3 - 8 {tt./ tt./ tt/) + 12 (Dtt) (tt^tt^) + 2 (Dtt) (7t ;ai tt") i? 

+ 4 (tt ;/1 TP"" 7r.„) i? + 8 (Dtt) (t^ ^ 7r ; „) - 4 (tt ; a ^ ;A ) (vr ;p , i?^) 

- 16 (tt, p TP"" R vp tt") - 8 (vr ;p 7r. v tt, p<j R ppua ) . 

where R pv and R pva p are the Ricci and Riemann tensors correspondingly. 

In what follows it will be useful to introduce dimensionless quantities in order to simplify 
formulae. Having in mind that we will apply our formalism to the problem of accretion, a convenient 
rescaling is 

x M -> r g x p , tt -)• Cr g TT, (5) 

where r g = 2M is the gravitational radius of a black hole, and the constant C (with mass dimension 
2) will be associated with the cosmological value of dtir, see below. 



III. CASE£ 2 +£ 3 

In this section we consider the accretion of a galileon, having nonzero £2 and £3, while other 
terms set to zero. This type of action (up to coefficients in front of £2 and £3) appears as the 
effective scalar field action, left over after taking the decoupling limit in DGP model of gravity. In 
the rescaled units the action reads, 

(6) 



2 In 0] the energy-momentum tensor was denned as T M ' y = l/s/—g {5S / Sg^v). 



r g C< 



d^Xy 



{dirf + k {dirf Utt 
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where e = 0, ±1, k = Cc^/r g and we let k to be positive or negative. Positive e corresponds to the 
canonical kinetic term, while positive both e and k correspond to the decoupling limit of DGpH. 
Note that the the constraints from [l[ imply e > and K > in the asymptotically flat space-time. 
In this case, for the static solution far from the source, tt' oc r~ 2 (such that tt' > 0), and the 
perturbations 5tt on the background tt propagate superluminally. In this paper, we, however, do 
not restrict ourselves to only positive e and n. 
The equations of motion derived from ([6]) are 

= 0, j M = 2e vr iM + n (ztt^Utt - d„ (c^) 2 ) , (7) 



or, 



e Dtt + k ( {air) 2 - (VVtt) 2 - R^tt^tt^ ) = 0. (8) 



We will also need the equation for perturbations 6tt on a non-trivial background n(t,x), in the 
limit of high frequencie^. From ([8]), we find the equation, 

G^V^VJtt = 0, (9) 

where 

GP" = (e + 2«D7r) g" u - 2kV^V u tt. (10) 

Note that in deriving ([9]), ()10p we neglected the term ~ dirdirSR, although there is a contribution 
of this term in the effective metric (|10p . see 0], However, this term is suppressed by Mm 0], so 
we can safely neglect it in our calculations (see also the discussion in Sec. 5). 

The propagation vectors for small perturbations can be found from the relation (see, e.g. [le 

<WY = o, (ii) 

where G^ v is the inverse matrix to G^ v , G^ V G^ V = 1. 

Since we are interested in solutions of a scalar field including some region inside the Schwarzschild 
horizon, we will be needing a regular (on the Schwarzschild horizon) coordinate system. E.g. we 
can take ingoing Eddington-Finkelstein (EF) coordinates. The EF coordinates (v, r) are connected 
to the usual Schwarzschild coordinates (t, r) in the following way, 



v = t + 

where 



dr 
J" 



/ = !--, (12) 
r 

in the rescaled units. The metric of a Schwarzschild black hole metric in EF coordinates reads, 

ds 2 = -fdv 2 + 2dvdr + r 2 dQ, (13) 
To describe a steady-state accretion we take the ansatz, 

/dv 
— +i>(r). (14) 



3 Note that the overall sign in the action is opposite to [J], so that in our notations the canonical kinetic (no-ghost) 
term is +(0ty) 2 . 

4 Not to be confused with "perturbation" which sometimes used in a different context: as a deviation of the solution 
from the homogeneous configuration due do the presence of a static source. 
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Note that by taking the ansatz (|14p we used the freedom in choosing the rescaling in ([5]), thus we 
took C to be equal dtir at the spatial infinity, 

C — dt7r\ r — OQ — v 7r\ r= 

thus setting the coefficient in front of v in (fTil) to unity. The term / f^dr in (JTIj) was introduced 
for convenience, such that for the homogeneous solution (with no black hole in the Universe), 
ip{r) = const. Since the current depends only on r, Eq. ([7|) can be integrated once to give 

r 2 f = A, (15) 

where A is a constant determining the total flux (to be fixed below). For the ansatz (114p the 
r-components of the current reads, 

f = 2eftf + k (-£ + //V 2 + . (16) 

Eqs. (|15p . (I16p can be obtained in a different way, namely, by integrating the equation for the 
energy-momentum conservation, T^ u = 0, which gives r 2 T£ = const, where for the ansatz (|14p 
one can find T£ = j r . The Eqs. (|15p and (|16p form an algebraic equation on tp'. The solution 
contains constant A as a free parameter, ip' = ip'(A,r). The physical solution will be chosen by 
the requirement that it is neither singular at the Schwarzschild horizon nor at the sound horizon. 
This seemingly "weak" requirement choses the solution uniquely, which in hydrodynamics is called 
transonic branch. To fix the transonic solution one needs to find a position of the sound horizon, 
i.e. one needs to study perturbations on the background solution 7r^(w,r). There are generically 
two solutions of (fl"5j) and (fTBT) . 

e r 2 / ± Ve 2 r 4 / 2 + nr (Af + jjg/Q {rf + 4/) 
^ ~ «r/(r/' + 4/) ' ( 17 ^ 

where the subscript (2, 3) means that the solution was obtained for the theory with £2 and £3 
terms in the Lagrangian. Below we consider solutions different cases corresponding to the choice 
of coefficients in the action ©. 



A. Canonical kinetic term, e = 1, K = 0. 

The standard case is recovered from (|1T[) by setting e = 1 and taking the limit k — > for the 
solution with the minus sign (alternatively, one can solve (|15p and (|16p setting k = 0), 

tinW = 4*]- ( 18 ) 



It is not hard to see that the physical solution corresponds to the choice A = 2 in (|18p . The sound 
horizon r* coincides with the Schwarzschild horizon, r* = 1, which is not surprising, since the speed 
of excitations for the the canonical term is the speed of light. The solution for n'(r) is shown in 
FigEl 



B. No canonical term, e = 0, k^O. 



It is interesting to look at the case when the quadratic kinetic term is absent. This case is 
simpler than the generic case (which we consider below), and in meantime it gives an insight to 
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the full problem. From Eq. (|17p . setting e = and substituting f(r) from (|12|) . we obtain, 



,/ , is 1 r\A(r - 1) + Kr] , . 

= 1 (4 ,._' 3 ) • (19) 

For the solution with the plus sign the singularity coming from the denominator of (|19j) at r = 1 is 
cancelled by the term (— f f~ l dr} in (|14p for any A; while the solution with the minus sign leads 
to a divergency of physical quantities. Another potentially dangerous point is at r = 3/4. The 
only way to avoid the singular behavior at this point is to to choose A = 3k. Thus, 

V>(„ocan) = (20) 

is the physical solution, in Fig. Q] the corresponding solution for ir' is shown. Let us now analyze 
propagation of perturbations. From (|10p setting e = one can find, 



O oo _ _ 2K ( / + /) , go: _ „ ( JV _L + *Jfy a n = ^M±*D£. (21) 

Then solving the equation (jlip with the covariant metric G^, being the inverse of (|2ip. we find 
two characteristics of the differential equation, rj = dv/dr, 

(2 ± \/3) r 

^ 3) " (i + ^) (V^rTs? 

The divergence of one the characteristic 77^3) (with the minus in the denominator) at r = 3/4, 
indicates the position of the sound horizon, r* = 3/4. Note that in this case the sound horizon is 
inside the Schwarzschild one. The other propagation vector, which corresponds to the perturbations 
being sent into the black hole, is finite everywhere at < r < 00. 



C. Generic case, e 7^ 0, k 7^ 0. 

Let us now turn to the generic case, i.e. when both terms in ([6]) are nonzero. First we consider 
the quadratic term in action (|6|) to be canonical, e = 1. The solution with the plus sign ()17|) . for 
any A, at infinity is approximated by 

v4 3 )+ = -^ + 0(i), 

thus violating our assumption on the homogeneity of the solution at the spatial infinity. Meantime 
it is regular at r = 3/4. The other solution, on the contrary, is well-behaved at infinity, 

V/ (2 ,3)_ = ^ (A + K ) + i (I + k) + O (r- 4 ) . (22) 

However, the solution (|22p diverges at r = 3/4 along with the characteristics ?7(2,3) f° r the generic 
choice of A. Thus, similar to the standard problem of fluid accretion, one needs to find such 
a value of A that at some point r = r* the solutions V'(2,3)± an d their derivatives match. The 
surface r = r* (called the critical point), is analogous to the acoustic horizon, from inside of which 
perturbations of tt cannot escape to the asymptotically flat regions (see below). At the critical 
point one of the propagation vectors 7/(2,3) diverges, i.e. in our case G vv = 0, implying G rr = 0. 
However, up to the non-relevant conformal factor, G rr is also being the coefficient in front of tp" in 
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0.0 rr 




FIG. 1. In the case C ~ £2 + £3, solutions for 7r' = — 1// + ^'( r ) are shown for different parameters of the 
model. The positions of sound horizon are shown by the dots. In case k < the sound horizon is outside 
the Schwarzschild horizon, while for k > the sound horizon is inside the Schwarzschild horizon. 



the equation of motion for the background solution. In order for the solution to be smooth at the 
sound horizon, the rest of equation of motion for ip (containing ip, ip' and terms) should also be 
zero at the sound horizon. Thus the procedure is as follows. Differentiating ([15)) once, we obtain 
the second derivative equation on ip. Then one can express ip" in terms of ip' (and there is no A in 
the expression, because it disappears during the differentiation). Now, the critical point is where 
both denominator and numerator are equal to zero. Solving simultaneously these two equations, 
one obtains an algebraical equation P(r) = 0, where P(r) is a polynomial 6th order in our case. 
The position of the sound horizon, r*, can be then found numerically from this equation. After 
finding r*, the constant of integration A is fixed. There are generically two possibilities, one of 
them corresponds to (physical) inflow and the other to (unphysical) outflow. 

For a positive k, corresponding to the decoupling limit of the DGP model, the sound horizon is 
inside the Schwarzschild horizon. The latter could be anticipated, since in this case the perturba- 
tions of the scalar field on nontrivial background of ir are superluminal. There is a regular solution 
in the whole range of r, < r < 00, see Fig. [TJ For k < 0, on the contrary, the sound horizon is 
outside the Schwarzschild horizon. The solution is smooth outside of the Schwarzschild horizon, 
and it becomes singular at some point r < 3/4, which can be seen by analyzing the square root in 
(|17p. This solution is nevertheless physical, since the pathological behavior happens inside of both 
acoustic and Schwarzschild horizons. 

After fixing the critical inflow solution, one can analyze the characteristics of the equations for 
perturbations, ^(2,3) and to check that r* is indeed the sound horizon. Fig. Q] shows behavior of the 
solutions for ir' for various cases. 

Note that one can also relax the requirement of positivity of e (being the requirement that the 
Lagrangian is having healthy canonical term) and to consider the case when the canonical kinetic 
term has a negative sign (ghost-like). Such a Lagrangian was, e.g. considered in the cosmological 
context in 0]. In our study, this situation can be easily reduced to one of the studied cases. Indeed, 
if both C2 and C3 are negative, then action is equivalent to the case C2 > 1, n > 0, up to the total 
sign in front of the action, so that the equations of motion do not change. The action with C2 < 0, 
C3 < 0, is equivalent to the studied above case C2 > 0, C3 < 0, up to the total sign in front of 
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the action. Note that we do not address the problem of ghosts here, we are only interested in the 
existence and stability of the classical solutions. 



IV. CASE OF NON-ZERO £ 4 



In this section we consider the case of nonzero £4, 



= AC 2 I d^x^-g 



, mf + k (9vr) 2 ( 2 (Ott) 2 - 2 (7r ;AW ,7r ; H - \v»R 



where k = c^C 2 /r 2 . In this case the problem of accretion becomes quite cumbersome. To start 
with, let us first consider the case when only £4 term is present in the action ([!]). 



A. No canonical term, e = 0. 

One can use the fact that the energy-momentum tensor for £4 is conserved, T^ u =00]. Since 
Tv does not depend on v explicitly for the ansatz (fT4"j) . we obtain, 

% = % (23) 
where from (|4"cj) . applying the ansatz Q14|) . after lengthy but straightforward calculations we find, 

T v = -fr (-rWV 2 + r 2 f" - 6r/' (/V 2 - l) - 4/V 2 ) , (24) 

where without loss of generality we have set the coefficient in front of £4 to unity. Now, the 
problem of accretion for non-zero £4 (and £« = for i ^ 4) is reduced to the solving (f23l) and (licj) 
and choosing the physical solution, if any. Combining (|23p and (|24p we obtain three branches of 
solutions ip'uy For A < or A > 8 there is no well-behaved solution. It is interesting that for any 
A in the range < A < 8 a smooth solution exists for ir' , covering the whole space, < r < 00, 
with corresponding characteristics finite everywhere, see Fig. [2j The critical point is absent and 
the solution at infinity is, 

A 1 ' 3 2 + A 2 ' 3 „/l\ 

The critical point exists for A = 8 at r* = 1/3; one can check that the sound horizon is indeed at 
r = 1/3. The solution critical solution is then, 

^4 



"(4),* - r _ J' 
which implies tt^ ^ = 0, see Fig. [2J 



B. Generic case, e / 0, k ^ 0. 



Let us now consider the case when besides the term £4 in the action, there is also the canonical 
kinetic term, £ ~ £2 + £4. Depending on the signs in front of £2 and £4, one can distinguish, in 
principle, four different cases. Let us first take £2 to be canonical, i.e. C2 = 1 and also C4 > 0. 
In this case a family of smooth solutions exists for a certain range of positive A, < A < A*{R). 



0.01 0.1 1 10 100 

r 

FIG. 2. Case, C <~ d- Solutions for it' = — 1// + ^'( r ) are shown for various values of A. The critical flow 
is for A = 8 (solid thick line), the critical point is indicated by the black dot, r* = 1/3. For the cases A = 4, 
1, the solution is shown by dotted lines, from top to bottom. 




FIG. 3. Case, C ~ £ 2 +C 4 - Solutions for tt' = -l/f + tp'(r) are shown for different R: R = -100, -1, -0.04 
(thick lines, from the top to bottom) and R — 0, i.e. only canonical term (thin line). The positions of sound 
horizon are shown by the dots, for R < they are outside the Schwarzschild horizon. For R < solutions 
are singular. 



For negative A the solutions diverge at the Schwarzschild horizon. The critical solution, realized 
for A = A*(k), is regular for < r < oo with the critical surface inside the Schwarzschild horizon. 
Interestingly, however, that even though smooth background solutions exist for < A < A*(k), the 
propagation vectors are nonexistent for some radii outside the Schwarzschild horizon. In fact, the 
equation of motion for perturbations become elliptic, indicating that solutions are catastrophically 
unstable. 

On the other hand, for e = 1 and R < 0, a critical solution exists, with well-behaved propagation 
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vectors, at least up to the sound horizon. The asymptotic form of the physical solution is, 



Note that the leading asymptotic behavior at infinity is the same as in the case of canonical term, 
up to the difference in constant A, which is simply the reflection of the fact that at infinity the 
term £4 becomes unimportant compared to £2. 



The results of the paper can be summarized as follows. For the Lagrangian Q, £ ~ £2 + £3, 
we have found, 

• case C2 = (canonical kinetic term is absent). Critical solution exists, Eq. (|20p . with the 
critical point at r* = 3/4, see Fig. [H upper (red) line. 

• case C3/C2 > 0. Critical solution exists, with asymptotic form (|22p at infinity, the critical 
radius r* is inside the Schwarzschild horizon. 

• case C3/C2 < 0. Critical solution exists as well, with the asymptotic form (|22p at infinity, 
however the critical radius r* is ouside the Schwarzschild horizon. 

For the Lagrangian £ ~ £2 + £4, 

• case C2 = (no kinetic term). There is a family of solutions parametrized by the value 
total influx, A. These solutions and the corresponding propagation vectors are smooth for 
< r < 00. Critical solution also exists for a particular value of the influx (for A = 8), with 
the critical point at r* = 1/3, see Fig. [2] 

• case C4/C2 > 0. A family of smooth solution exist for some range of positive A. The 
critical solution also exists with a critical point r* < 1. However, for the found solutions the 
propagation vectors do not exist in some region of r, r > 1. 

• case C4/C2 < 0. There is the well-behaved critical solution and propagation vectors, cor- 
responding to the solution. The sound horizon is outside the Schwarzschild horizon, see 



It is interesting to see if we get any new constraints on the Lagrangian compared to those ob- 
tained for the non-covariant version of galileon. Since in our study the space-time is asymptotically 
flat, it would be natural to apply constraints on a obtained in [l[ for flat space-tim^|. For the 
Lagrangian containing only Co&nd £3 we did not find any constraints on the coefficients C2, C3. 
Note that C2 > was found in [ill by requiring the model to be ghost-free. Since we did not address 
the issue of ghosts in our studjo, it is not surprising that we did not find any constraints on C2. 
It is more interesting, that if we impose positivity of C2 we still do not get any constraints on C3 
(meantime the study in [1] implies that C3 > if applied to our case). There is a subtle point 
here: for our solution the case C3/C2 < corresponds to the critical point outside the Schwarzschild 
horizon, and formally the solution diverges at some point inside both horizons. This behavior, 




V. RESULTS AND DISCUSSION 



Fig. El 



5 Then in [l|] one should make substitution di — > a. 

6 In fact, the presence of ghost does not affect the solutions we found, since the overall sign in front of the action is 
unimportant. 
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however, cannot be used to invalidate the solution, because the region inside of both horizons is 
causally disconnected from the asymptotically free region, so, in principle, inside the Schwarzschild 
horizon the solution can change to another one (the one which does not satisfy the steady-state 
ansatz). 

For the Lagrangian £ ~ £2 + £4 we also obtained an interesting result. In particular, for both 
C2 and C4 positive (and C3 = 0), the solution for the steady-state accretion is unstable, because 
the propagation vectors do not exist in some interval of r outside the Schwarzschild horizon. On 
the other hand, for C4/C2 < 0, we obtained a well-behaved critical solution with corresponding 
propagating vectors. These results are in contrast to those obtained in [l| for a static configuration 
of non-covariant galileon, where for C3 = and asymptotically flat space-time the coefficient C4 has 
to be positive or zero. 

The propagation of small perturbations for galileon is superluminal for a range of parame- 
ters. Thus, it is not surprising that in this range of parameters the sound horizon is inside the 
Schwarzschild horizon, implying that one can in principle communicate with the interior part of 
a black hole, r > r*. The same effect was found for "superluminal" k-essence Similarly 
to k-essence accretion, the acoustic space-time for accreting galileon is globally hyperbolic, in 
spite of the superluminal perturbations. This can be explicitly demonstrated by constructing a 
Schwarzschild-like metric of the perturbations of the field. We take the case £ oc £3 for simplicity. 
Substituting (|13p and (|20p into (|10p one can find the effective contravariant metric. By inverting 
the obtained metric, one finds the covariant effective metric G^ u . Making successive change of 
coordinates, dv = dT — Gt r /Gudr, r = i? 3 / 4 and T — > y/l2T, we finally get, 

dS 2 = -(4R 4 / 7 - 3)R 2 / 7 dT 2 + —£L-—dR 2 + R 2 dfl. (25) 

The above expression for the effective metric is similar to the Schwarzschild metric, with the 
important difference that the position of the event horizon is now at i? 4 / 7 = 3/4. 

Although the global hyperbolicity is not affected by the presence of the superluminal propa- 
gation in this case of steady-state accretion of galileon, it is still unclear, however, how to embed 
the existence of the new "event" horizon into the standard thermodynamics of black holes. In 
particular, since the sound horizon is inside the Schwarzschild one, the latter loses its meaning of 
the universal event horizon, therefore one may need to modify the entropy assigned to a black hole. 

It is interesting to note, that for a range of parameters, our solutions for accretion possesses 
the analogue of the Vainshtein behavior, similar to the case of static sources|f| Indeed, let us take 
the case £ ~ £2 + £3 for definiteness and e > 0. In the limit r — > 00 the solution is given by (|22|) . 
which has the same form as the solution for the standard kinetic term (|18p . apart from the overall 
constant A — > A + k. In this regime the kinetic term is dominant, and assuming that non-minimal 
coupling is given by a term ~ irT/Mp (in the Einstein frame) a test particle would experience an 
additional force, F n , which is to be compared to the gravitational force -F gra v; 

F n V7r r 2 •h 00 r g {A + K) 



(26) 

-Fgrav Mp T g Mp 

where we restored the physical units. On the other hand, if also k > 0, then the solution changes 
its behavior at the radius, 



r g [4k(A + k)\ 



1/3 



7 A similar arguments can be repeated for other models of galileon, considered above. 

8 Note, however, that in order to observe the Vainshtein effect, the matter must be coupled to the galileon non- 
minimally, so that the fifth force (due to the direct coupling) leads to the deviation from General relativity on 
large scales, while on small scales the galileon is screened. Meanwhile, our results for the accretion do not depend 
on coupling to matter, since the solution is matter-free. 
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where r* is a analogue of the Vainshtein radius. The "accretion" Vainshtein radius is related to 
the "usual" one, ry, defined for a static source as follows, 



\ V3 

r* / ^"oo^o 



Thus, for the range of distances r g < r < r*, the ratio of the galileon and gravitational forces is, 



F n Tr 00 r g (A + K) ( r 



F„ rav M P \r 



3/2 



which has an additional small factor (r/r*) 3 / 2 compared to (j26j) . being a manifestation of the 
Vainshtein effect. This can also be seen in Fig. (pQ): the behavior of ir' for k > is smoother closer 
to the horizon, compared to the canonical kinetic term. Note, that if e and k have different signs 
(e.g. e > and k < 0), then solutions do not exhibit the Vainshtein behavior, due to the fact that 
the term in ()17p . which responsible for this (the second term inside the square root) is negative 
and therefore it can never be dominant. A similar analysis can be made for the Lagrangians of the 
form £ ~ £2 + £4, and one can find that the solutions have the Vainshtein behavior in this case 
as well. 

It is worth to mention the change of the black hole mass during the accretion of the galileon. 
Since by definition the problem we studied was steady-state in the test-fluid approximation, the 
rate of the change can be found via the total flux at the infinity, r — > 00. In the Schwarzschild 
coordinate system the total flux is oc (r 2 T[). Expressing T[ through the components in the EF 
coordinates, we obtairjfl, 

irf"- 4 ^- (27) 

where the physical units have been restored. In particular, for the case when only the canonical 
term is present, A = 2, one gets the known result (compare, i.e. with [13], the additional factor 2 
in the above expression is due to the choice of the canonical term (dir) 2 , instead of l/2(<9"7r) 2 ). The 
flux can be easily made negative by changing of the overall sign of the Lagrangian for tt, so that the 
black hole will loose mass instead of gaining it. Usually, such a change of sign is accompanied by 
ghosts. However, e.g. in Q] it was shown that even though the £2 term is ghost-like, nevertheless 
the full Lagrangian, £ ~ £2 + £3, is ghost-free in the vicinity of the cosmological attractor. One 
can think of exploiting this particular feature to construct negative flow onto a black hole. In this 
case, however, one immediately stumbles the following problem. In our study the space-time is 
taken to be asymptotically flat, implying that the term £3 is subdominant in comparison to £2, 
implying the appearance of ghosts in the model studied in 0]. (Note that the expansion of the 
universe in the attractor is the key to to avoid ghosts in [9]). Thus, in the region of space-time 
where the gravitation of a black hole is stronger than the cosmological expansion, one may naively 
expect the appearance of ghost. This question, however, deserves a separate and more careful 
study. 

Let us in conclusion discuss the assumptions we made to find the above results. First of all we 
assumed that the accretion is steady-state, which is physically reasonable when the Lagrangian is 
shift-symmetric and it is enough time for the system to evolve before the steady-state is established. 
We also neglected the backreaction of the accreting fluid on the metric, since the backreaction can 
be made parametrically small. Indeed, the total flux onto the black hole is proportional to r 2/ k 2 <1 



The same expression can be obtained by calculating the flux at of the Killing energy across the Schwarzschild 
horizon, dM/dt — Area x T vv . 
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(f2?|) . In the limit C3 = const, r g — > 0, tt^ — > Ofand keeping ■n (yo lr g = const, so that the solution, 
expressed in dimensionless units is not affected 10 !) . the total flux goes as ~ r g , while the mass of the 
black hole ~ r g , thus validating our the assumption. Another source of possible inaccuracy in our 
calculations is that we also neglected the kinetic mixing term of perturbations and gravity. E.g. 
the term ~ dTrdtrSR was neglected in (jlOp . The contribution of this term is to compared to the 
contributions of the terms left over in (110p , one can check that in order to neglect the backreaction 
of the perturbations, one needs (ddir) 3> C3 (dir) 4 /Mp, which again can be maid parametrically 
small by taking the limit we used before. 

To summarize, in this paper we studied the steady-state accretion of the covariant galileon 
onto a Schwarzschild black hole. Subject to the assumptions we made, our study of existence and 
classical stability of solutions did not show any constraints on the parameters of the Lagrangian 
of the form C ~ £2 + £3. In the event when £4 is present, the constraints on the parameters of 
the theory we obtained are in a sense orthogonal to the constraints found in other studies. We 
also notice that for some range of parameters of the theory signals can be sent from the inside of 
the Schwarzschild horizon, being a consequence of superluminal propagation of perturbations. The 
accretion of the energy flow can be made negative by appropriate choice of the coefficients in the 
Lagrangian. 
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